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We report on a recent extraction of the higher twist contributions to the deep inelastic 
structure functions F^*'''' (x,Q'^) in the large x region. It is shown that the size of 
the extracted higher twist contributions is strongly correlated with the higher order 
corrections applied to the leading twist part. A gradual lowering of the higher twist 
contributions going from NLO to N'^LO is observed, where in the latter case only the 
■ ^ ■ leading large x terms were considered. 

II ■ In wide kinematic regions the deeply inelastic structure functions can be described by their 

C1h| leading twist contributions within Quantum Chromodynamics (QCD). Higher twist correc- 

tions [1] emerge both in the region of large [2,3] and small values [4,5] of the Bjorken variable 
X. The leading twist sector, both for unpolarized and polarized deeply inelastic scattering, 
is well explored within perturbative QCD up to the level of 3-loop, resp. 2-loop, correc- 
tions [6,7], including the heavy flavor contributions [8]. On the other hand, very little is 
known on the scaling violations of dynamical next-to-leading twist correlation functions and 
^1 I the associated Wilson coefficients [1], even on the leading order level. In many experimental 

C\J ■ and phenomenological analyzes, cf. [2,3], higher twist contributions are parameterized by 

^D \ an 'Ansatz' [2], which is fitted accordingly. Within QCD this ad-hoc treatment cannot be 

r~^ . justified, performing at the same time a higher order analysis for the leading twist terms. 

^^ ' Since neither the corresponding higher twist anomalous dimensions nor Wilson coefficients 

OO , were calculated, the data analysis has to be limited in the first place to the kinematic domain 

in which higher twist terms can be safely disregarded. 

In the following we report on the determination of the higher twist contributions in 
the deeply-inelastic structure functions Fj^'"^ {x,Q'^)., see Rcf. [9] for details. Higher twist 
contributions were also studied in deep- inelastic neutrino scattering, cf. [10]. Also in the 
C^ ' case of polarized deeply inelastic scattering higher twist corrections are present in general. 

Since the polarized structure functions are measured through an asymmetry, the effect of 
higher twist contributions in the denominator function has to be known in detail. In [11] 
no significant higher twist contributions were found. Other authors claim contributions in 
the low X region [12], which is also the region of very low values of Q^ . These effects need 
further study. 

In the case of the structure functions F!^^''^ (x, Q^) we investigate the flavor non-singlet 
contributions in the large x region for Q^ > 4GeV^,iy^ > 12.5 GeV^, cf. [13]. One gener- 
ally may consider flavor non-singlet combinations and perform a three-loop QCD analysis, 
which requires the 0{al) Wilson coefficients [15] and the 3-loop anomalous dimensions [6]. 
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^For power correction analyzes in the resonance region see [14], Here the concept of the twist-expansion 
is not applicable, except assuming duality. 
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The analysis can even be extended effectively to 4~loop order, since the dominant contribu- 
tion there is implied by the 3-loop Wilson coefficient [7], parameterizing the yet unknown 
4-loop anomalous dimension with a ± 100 % error added to an estimate of this quantity 
formed as Pade-approximation out of the lower order terms. A comparison with the 2nd 
moment of the 4-loop anomalous dimension calculated in [16] showed [13] that the agree- 
ment is better than 20 %, which underlines that the above approximation may be possible. 
We limit the QCD-analysis of the twist-2 contributions to this representation since neither 
as(^^) nor the splitting and coefficient functions are known beyond this level. Furthermore, 
heavy quark and target mass corrections are applied. 

The evolution equations are solved in Mellin-A^ space, cf. [17]. The non-singlet structure 
function at the starting scale of the evolution, Qq, is given by 

oo 
k=Q 

with C^^{N) the expansion terms of the non-singlet Wilson coefficient with Cq{N) = 1, 
a,s(Q^) = '^s(Q^)/(47r) and /f' ' {N,Q^) the corresponding combination of quark distri- 
butions, cf. [13]. Here we identify both the rcnormalization and factorization scale with 
Q^. Beyond 0{al) dominant large x contributions to the Wilson coefficient were calculated 
in [18]. 

We then extrapolate the results to the region 4 GeV^ < W^ < 12.5 GcV^ and determine 
effective higher twist coefficients Cht{x, Q^) given by 



Fr{x,Q^) = Fr\x,Q^) 



OTMc[Fr^{x,Q^)\ Cht(^,Q^) 



(2) 



Here Otmc[ ] denotes the operator of target mass corrections. 

QCD corrections beyond N'^LO are known in form of the dominant large- a; contributions 
to the QCD-Wilson coefficients [18]. Since these corrections do quantitatively only apply 
in the range of large x we do not use them in the twist-2 QCD-fit, because here the data 
are mainly situated at lower values of x where beyond 4--loop order other contributions to 
the Wilson coefficients, which are not calculated yet, are as important. Furthermore, the 
4-loop anomalous dimensions are yet unknown beyond the 2nd moment. The leading large 
X contributions are given in terms of harmonic sums of the type 5'i^i_...^i(iV) which obey a 
determinant representation [19] in single harmonic sums Si{N). The effective higher twist 
distribution functions C^'rp(x) extracted are shown in Figures 1 from NLO to N^LO. Here 
we averaged over the values in Q^ within the x-bins. The leading twist terms are those 
given in [13], with the values of A^^^^^ = 265 ± 27, 226 ± 25, 234 ± 26 McV, resp. in NLO, 
NNLO, and N^LO. Both for the proton and dcuteron data Cwt{x) grows towards large 
values of x, and takes values ~ 1 around x = 0.6. The inclusion of higher order corrections 
reduces Cwt{x) to lower values with a gradually smaller difference order by order. Yet for 
the highest bins, x > 0.8, the effect of the large x resummation terms is important. Earlier 
higher twist analyzes [3] limited to the next-to-leading order corrections are thus corrected 
by factors of 2 and larger at large x to lower values. Soft resummation beyond NLO was 
considered in [20]. In the present analysis we limited the investigation to the inclusion of 
the large x terms in N^LO which arc still in the vicinity of a nearly complete QCD analysis 
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as outlined above. The present description is likely to be final for values of a: < 0.8. Beyond 
this range there are only few data. More data in this interesting region would be welcome 
and can be obtained at planned high-luminosity colliders such as EIC [21]. 
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Figure 1 : The coefficient CHT(a^) for proton and deuteron data in the large The curves cor- 
respond to the cases of twist-2 corrections in NLO: dotted line, NNLO: dashed line, N'^LO 
dash-dotted line, and asymptotic N'^^LO: full line, cf. [9]. 

In the present analysis we extracted the large x dynamical higher twist contributions to 
the structure functions in a model-independent way. It would be interesting to compare 
moments of the term Cwvix) to lattice results, which allow to simulate the moments of the 
corresponding higher twist correlation functions, in the future. 
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